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Using Kohn-Sham wave functions and their energy levels obtained by density-functional-theory 
total-energy calculations, the electronic structure of the two c(2 x 2) phases of Na on Al (001) 
are analysed; namely, the metastable hollow-site structure formed when adsorption takes place at 
low temperature, and the stable substitutional structure appearing when the substrate is heated 
thereafter above ca. 180K or when adsorption takes place at room temperature from the beginning. 
The experimentally obtained two-dimensional band structures of the surface states or resonances 
are well reproduced by the calculations. With the help of charge density maps it is found that in 
both phases, two pronounced bands appear as the result of a characteristic coupling between the 
valence-state band of a free c(2 x 2)-Na monolayer and the surface-state/resonance band of the Al 
surfaces; that is, the clean (001) surface for the metastable phase and the unstable, reconstructed 
"vacancy" structure for the stable phase. The higher-lying band, being Na-derived, remains metallic 
for the unstable phase, whereas it lies completely above the Fermi level for the stable phase, leading 
to the formation of a surface-state/resonance band-structure resembling the bulk band-structure of 
an ionic crystal. 

PACSnumbers: 82.65.M, 79.60.B, 68.35.B, 73.20.A 



I. INTRODUCTION 

The adsorption of alkali metal atoms on metal surfaces 
has attracted much attention in recent years partly due 
to the discovery of a variety of new adsorbate phases, 
in particular, structures that involve a reconstruction of 
the metal surface induced by the alkali metal atoms (see, 
for example, Refs.EJta and references therein). A common 
feature of these systems is that often there is a metastable 
phase at low temperature involving no reconstruction of 
the metal surface while at higher temperatures, the stable 
reconstructed phase occurs. 

In the present paper we present a combined theo- 
retical and experimental investigation of one such sys- 
tem, namely that of c(2 x 2) phases of Na on Al(OOl). 
Early dynamical-theory analyses of low energy electron 
diffraction (LEED) intensitiesclu concluded that Na occu- 
pied the four-fold hollow site. It was first demonstrated, 
by high resolution core-level spectroscopy (HRCLSp 
and surface extended x-ray absorption fine structure 
(SEXAFS)EI studies that the hollow site is taken only for 
preparations at low temperature (LT), and a different, 
stable structure is formed by adsorption at room tem- 
perature (RT) or by heating of the LT phase above ca. 
160K. It was shown later by a density functional the- 
ory (DFT) studyu and by a new LEED analysis!! that in 
the RT phase the Na atoms occupy substitutional sites, 
where every second Al atom in the top layer is displaced 
and a Na atom adsorbed in its place. The result of an 
x-ray photoelectron diffraction (XPD) studyd of the RT 
phase concluding that the surface contains two domains, 
could not be confirmed!! 

In what follows, we perform a theoretical analysis of 



the electronic structure of the two c(2 x 2) structures 
and compare the results with those of angle-resolved ul- 
traviolet spectroscopy^. The measurements have been 
performed in what we call "polar scan" modes which 
deliver displays of photoemission intensities as a func- 
tion of energy and wave vector component parallel to 
the surface, lying in selected symmetry directions (T-M 
and T-X) in the two-dimensional (2D) Brillouin zone (see 
Fig. 1). The displays yield directly 2D band structures 
of surface states/resonances (i.e. surface states or reso- 
nances, depending on the position of the states in or out 
of the gap of the 2D projected bulk bands), which may 
be compared with calculated results. 

The basis of our theoretical analysis are Kohn-Sham 
wave functions and their -energy levels obtained by DFT 
total energy calculations^. 2D band structures are de- 
rived and compared with the experimental results. The 
obtained, satisfactory agreement between theory and ex- 
periment may be regarded as a confirmation of the pro- 
posed atomic structure models mentioned above, and a 
useful basis for further studies of the properties of these 
surfaces. 

Also single-state charge density distributions of occu- 
pied and unoccupied states are derived from the DFT 
calculation, and used for analysing the character of the 
bands. We find that adsorption leads, in both LT and 
RT cases, to two main bands of surface state/resonances 
as the result of a coupling between the valence-state 
band of a free c(2 x 2)-Na monolayer and the surface- 
state/resonance band of the Al surface. Across the Fermi 
level the Al-derived band is shifted down, and the Na- 
derived band is shifted up. As a consequence, a charge 
transfer is taking place from the adsorbate directly into 
the pre-existent surface state/resonance of the substrate. 
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FIG. 1. Surface Brillouin zone of (1 x 1) (full lines) and 
c(2 x 2) (broken lines) structures. The symmetry points T, 
M, and X are referred to the (lxl) structure. 

For the LT (hollow) structure, it is found that the 
coupling is relatively weak and the Al-derived state al- 
most retains the perfect (lxl) periodicity of the clean 
surface. For the RT (substitutional) structure, the Al- 
derived state exhibits a clear c(2 x 2) periodicity of the 
reconstructed Al surface. The Na-derived band of the 
RT structure is completely empty and leads, except for 
the still existent background bulk-continuum, to a sur- 
face electronic structure having a character of an ionic 
monolayer lying on the surface. We discuss possible con- 
sequences of this finding. 

The paper is organised as follows: In Section II 
the experimental method is outlined and is followed in 
Section III by a description of the calculation meth- 
ods. Section IV contains the results for the surface 
state/resonance band structures and the analysis of the 
character of the bands in terms of single-state charge den- 
sity distributions. Also the charge transfer processes are 
analysed by using density differences. Section V contains 
the discussion and Section VI, the conclusion. 

II. EXPERIMENT 

The photoemission experiments!^ were performed in 
a VG ESCALAB Mark II spectrometer modified for 
motorized sequential angle-scanning data acquisitionliil, 
and with a working pressure in the lower 10 -11 mbar 
region. Ultraviolet photoelectron spectroscopy (UPS) 
measurements were done using unmonochromatized He I 
(21.2 eV) radiation from a discharge lamp. The 150-mm 
radius hemispherical analyzer was run with an energy res- 
olution of 50 meV. Contamination free surfaces were pre- 
pared by a combination of Ar + sputtering and annealing 
at 500 C. Na was evaporated from a carefully outgassed 
SAES getter source. Particular care was taken to en- 
sure ultraclean Na deposits: All parts of the evaporation 
source, except the tiny exit slit, were surrounded with 
liquid-nitrogen cooled walls. In this way, the pressure 
during evaporation could be kept as low as 2.5xl0~ n 
mbar. Na coverages were determined accurately (±0.03 
ML) from core-level photoelectron intensities^!. The 
sample temperature was measured with a thermocouple 
in mechanical contact with the sample holder. The tem- 



perature gradient between the sample and the sample 
holder was determined in a separate calibration experi- 
ment with an additional thermocouple spot-welded onto 
a dummy sample. Sample temperatures given here are 
corrected for this temperature difference and are esti- 
mated to be correct within ±10 K. 

In order to obtain two-dimensional band structures 
of surface states/resonances, polar scans along the T- 
M and T-X directions of the Al(001) surface Brillouin 
zone were performed, recording at each angular setting 
(A© = 2°) the entire photoelectron spectrum between 
—0.4 and 4.3 eV binding energy. The experimental data 
sets l(Ei,Q) acquired in this way were mapped onto a 
regular (J5j,k||) grid and visualized as gray-scale plots 
with low intensities in white and high intensities in black. 

III. CALCULATION METHOD 

The ab initio DFT total energy calculations, agd com- 
parison with LEED results are described in RcfsfrH. The 
calculations were performed using the local-density ap- 
proximation (LDA) for the exchange-correlation func- 
tional. jEurther details about the method can be found 
in RefsfltM 

The following two kinds of diagrams have been used in 
the present analysis: 

1. 2D band structures 

For deriving 2D band structures we use the follow- 
ing procedure: 

The projected density of states (DOS) onto a cho- 
sen atomic orbital, ?pAo( r ), is defined as 

JV AO (e) = X) JV Ao(k||,E), (1) 
k n 

where the "k|| -resolved" projected DOS is given by 
iV AO (k||, £ )=^| [ d 3 r^ (r^ k|1 ,e'(r)r 



7T (e' — e) 2 + 7 2 

We have used for ^ao ( r ) the eigenfunctions of the 
isolated pseudo-atoms from which the pseudopo- 
tentials were derived for the total energy calcula- 
tionUl. The integral was truncated at a cutoff ra- 
dius r c . (We used here r c = 3.7A .) In the standard 
supercell method using a slab geometry, a Kohn- 
Sham state, Vk|i,e( r )j can be specified, except for 
the presence of degeneracy, by indices k|| and e 
(both discrete), the parallel wave vector and the 
energy, respectively. The lifetime broadening pa- 
rameter, 7, is introduced for convenience in numer- 
ical work and for improving visibility of peaks in 
the resulting plots. (We take here 7=0.5 eV.) 
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A simple sum of iVAo(k||, e)'s are f° rme d over 
atomic orbitals of a specified atom. For example, 
for Na we have, 

A Na (k||,e) = iVNa3 S (k||,e) + A Na 3px (k|| , e) 

+ ^Na3p H (k||,e)+7VNa3 P ,(k||,e) . (3) 

The corresponding quantities, iVAi(kii, e), are eval- 
uated for the Al atoms in the uppermost layer. 
The maximum peak of these quantities was always 
found at the bottom of a surface state/resonance 
band lying at T. Those peaks having a fractional 
ratio to the maximum higher than a properly cho- 
sen fixed value ("Min. Fraction" in Tab. I), are 
selected and their positions in (kii , e)- space dis- 
played, using squares for Na and circles for Al (see 
Tab. I). The plots along X-T and f-M (see Fig. 1) 
are combined, and are presented in Figs. 2 and 3. 
We see that they produce satisfactorily 2D band 
structures to be compared to the experimentally 
obtained ones. 

The bulk-band continuum is separated into discrete 
bands due to the use of finite (nine layer) Al slabs. 
They appear in Figs. 2 and 3 as weak features. The 
discreteness of bulk bands may have caused a small 
energy shift of the surface resonance bands lying in- 
side the bulk continuum, because these are replaced 
by one of the discrete bulk bands lying nearest to 
them. 

We note that in Figs. 2 and 3 the theoretical bands 
for the LT and RT structures repeat themselves in- 
side the (lxl) Brillouin zone, exhibiting the c(2 x 2) 
periodicity. That is, the range T-M (M being re- 
ferred to (1 x 1)) is halved, and the bands in the 
second half become a mirror image of theJirst half. 
We call it here "back-folding" (cf. RefEJ). It is 
a consequence of our use of density of states pro- 
jected to one specified atom in each c(2 x 2) surface 
unit cell. We note that this is obviously a theoreti- 
cal construct and not quite adequate for represent- 
ing fully the character of the wave functions of the 
bands, particularly in the LT structure. In fact, 
there is no information included about the relation 
between the values of wave functions around the 
two Al atoms in a unit cell (see Fig. 4). Thus, 
in reality the wave functions may happen to have 
approximately the (lxl) Bloch-type periodicity - 
the back-folding still occurs. The same would also 
result even for the clean-surface band if the arti- 
ficial c(2 x 2) unit cell would be imposed. In our 
back-folded band structure, a vanishing deviation 
from the (lxl) Bloch-type periodicity would be- 
come visible only in a vanishing band gap at a new 
Brillouin zone boundary, that is, in our case on the 
line halving the range T-M. 

On the other hand, the photocmission intensity 
would be determined by a matrix element 



M(k| h e)= / d 3 rV> f *(r)VV>k ll>E '(r). (4) 

The symmetry of the structure leads to the "se- 
lection rule" , which selects out initial-state bands 
according to the symmetry relation between the 
initial- and final-state wave functions. In the case of 
the LT structure, the final state ip{( v ) samples the 
initial state ipk n ,e'( r ) at both of the two Al atoms 
in a unit cell. The selection rule can take place dif- 
ferently for the first and second halves of the range 
T-M, destroying the mirror symmetry between the 
two halves. In particular, if the wave functions have 
nearly the (lxl) Bloch-type periodicity, the selec- 
tion rule leads to "unfolding" of the back- folded 
band structure. 

2. Charge density distributions 

The charge density is derived from Kohn-Sham 
wave functions as 

pw = EEk,^)i 2 . ( 5 ) 

ky eocc 

where e occ labels the occupied levels. We use p(r) 
for the study of charge transfer ocurring at adsorp- 
tion. 

The single-state charge density is defined here as 
Pk| h£ (r) = |V>k lh e(r)| 2 . (6) 

and used for analysing the character of surface 
states /resonances. 

We should keep in mind that the charge density dis- 
tributions are derived from pseudo-wave-functions 
resulting from the use of pseudopotentials, so that 
they give a correct distribution only outside the 
critical radius of the pseudopotentials. 

We use cross-sections of the charge densities on 
planes perpendicular and parallel to the surface. 
The perpendicular cross-sections taken are speci- 
fied in Fig. 4. 

IV. RESULTS 

In the next three Sections (A, B and C) we discuss 
successively the results for the clean Al (001) and the 
two c(2 x 2)-Na/Al (001) surfaces. It is to be noted that 
the Figures are arranged, independently of the text, in 
the way to facilitate a visual comparison between the 
different structures. 
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FIG. 2. Comparison of experimental (left panel) and calculated (right panel) 2D band structures. The symbols used for the 
calculated bands are specified in Tab. I. (a) Clean surface, (b) LT phase (hollow site), and (c) RT phase (substitutional site). 



A. Clean Al (001) 

Figure 2a shows the experimental (left panel) and theo- 
retical (right panel) band structures for the clean Al (001) 
surface. It can be seen that both exhibit a surface-state 
band which has a free-electron-like form. The paramaters 
obtained by fitting by parabolas, as indicated by broken 
lines in Figure 2a, are shown in Table II. The energy po- 
sition of the band at T, Eq, agrees well between theory 
and experiment, with values of 2.68 and 2.76 eV below 
the Fermi level, respectively. The theoretical pJpe lies in 
the range of the results of other calculationsESS: 2.6 - 
2.9 eV. The experimental value is in very good agreement 
with earlier studicsEjHj andEj. 



In Fig. 2 (left panel) it can be seen that the experimen- 
tal results display a number of additional features. Those 
which are common to all structures studied on this sur- 
face are assigned to bulk bands, as indicated by arrows. 
These bulk features do not appear in the calculated band 
structures (right panel), reflecting the fact that the den- 
sity of states projected onto the uppermost Al atoms are 
dominated by the surface states. 

In Fig. 6a we display the single-state charge density 
distributions at V for the main band. The two cross- 
sections (100) and (110) are defined in Fig. 4a. An im- 
portant characteristic of the surface state of the clean 
Al (001) surface is that the charge density shows a pro- 
nounced maximum just on top of the surface Al atoms. 
This can also be seen in the cross-section parallel to the 
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surface through the electron density maxima as shown 
in Fig. 7a. This particular feature of the surface state 
is apparently crucial in the formation of the electronic 
structure of the LT phase, as we see below. 
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B. LT c(2 x 2) hollow structure 

1. Al-derived band 

In Fig. 2b we show the experimental (left panel) and 
calculated (right panel) surface state/resonances for the 
LT c(2 x 2) hollow structure. The lower lying main band 
is Al-derived as indicated by the circles in the theoretical 
plot (cf . Table I) . The position of the calculated band at 
r lies somewhat lower in energy than the experimental 
value; compare 3.61 to 3.12 eV (see Table II). (Notably 
the calculated results of the earlier theoretical work of 
Chulkov and SilkintSis also «3.6 eV, although their non- 
self-consistent value of vertical distance 2.05 A, from Na 
to Al is different from ours.2.35 A, which is nearer to the 
experimental result 2.57 A0.) Compared to the position 
of the surface-state band of the clean Al (001) surface, 
the Al-derived band lies lower in energy, by 0.93 and 
0.36 eV, as obtained by the calculations and as deter- 
mined from experiment, respectively. The nature of this 
downward shift will be discussed below. The differences 
between theory and experiment may have been caused by 
the approximations used in the total energy calcttLation 
(use of LDA, Rydberg cut, k-point sampling etcJlMfJ). 

The experimental results show clearly that the Al- 
derived state does not have the c(2 x 2) periodicity, but 
rather the (1 x 1) periodicity of the clean surface. In fact, 
we find it a significant experimental result that the main 
band in the LT phase has, throughout the whole range of 
T-M, almost the same form as that of the clean surface, 
being only shifted down. 




FIG. 4. Geometries of the c(2 x 2)-Na/Al (001) surface 
structures: (a) clean surface (b) LT phase (hollow site), (c) 
vacancy structure and (d) RT phase (substitutional site). The 
positions of perpendicular cross-sections used for charge den- 
sity distributions are indicated. The white circles represent 
Na atoms and the grey circles Al atoms. 



FIG. 3. 2D band structures of (a) the LT phase (hollow 
site), (b) the RT phase (substitutional site), and (c) the va- 
cancy structure, including the energy range above Ef. The 
symbols used are specified in Tab. I. 



In Fig. 6b we display the single-state charge density 
distribution at T for the Al-derived band of the LT phase. 
The two sets of cross-sections differently chosen for Al 
and Na (cf. Fig. 6c) are defined in Fig. 4b. It can be 
noticed in Fig. 6b that the maxima of charge density re- 
mains on top of the surface Al atoms, almost unchanged 



5 



from those of the clean surface shown in Fig. 6a.— (This 
has also been pointed out by Chulkov and Silkinlia.) In 
the cross-section parallel to the surface shown in Fig. 7b 
only a small deviation from the (1 x 1) structure is seen, 
having a c(2 x 2) period, with some indication of a charac- 
ter of a bonding state between the Na and Al atoms. The 
smallness of the deviation explains the similarity of the 
experimental band structures between the clean surface 
and the LT phase. It is to be noted that in the pho- 
toemission processes from this band, not only the initial 
states, but also the final states have apparently main- 
tained approximately the same periodicity of the clean 
surface. The theoretical curves in Fig. 2b appear differ- 
ent, but this is only caused by the back- folding as already 
mentioned. 



2. N a- derived band 

In Fig. 2b dark spots can be seen at T and M, near the 
Fermi energy Ep, in the experimental results. Although 
from the figure presented here it is perhaps difficult to 
distinguish, we find in the experimental data that the 
dark feature at M for the LT phase is markedly different 
in character from that of the clean surface, in that the 
intensity for the LT phase is stronger and some tailing 
with dispersion is exhibited, as is the case also for the 
feature at T. From the calculated bands of the LT phase 
(right panel of Fig. 2b), we can see that there is a surface 
state/resonance band with an energy of about 0.7 eV 
below the Fermi level at T. As indicated by the open 
squares, it is a Na-derived band. We assign this band as 
giving rise to the experimentally measured intensities at 
r and M at ~0.4 eV below the Fermi level. The existence 
of this, band was first theoretically predicted by Benesh 
et altj and was also reproduced by Chulkov and Silkinta. 
In Fig. 3a we show the same calculated band structure 
as in Fig. 2b, but where the energy region extends higher 
into the positive range. Here we can see that the part 
of the Na-derived band above Ep exhibits marked band- 
structure effects in the middle of T-M, due to the c(2 x 2) 
periodicity of the adsorbed Na layer. 

In Fig. 5a we show the band structure of a free c(2 x 2)- 
Na monolayer for comparison. It can be observed that 
at r the lowest band lies «0.8 eV lower in energy as 
compared to the Na-derived band of the LT phase. Cor- 
respondingly the occupied part of the band is notably 
larger for the free monolayer. 

Figure 6c displays the distribution of the Na-derived 
surface state / resonance at T which clearly shows an anti- 
bonding character, the nodal line (not shown) running 
between the Na- and Al-layers. A comparison with 
Fig. 5b shows that the strongly smeared out character 
of the density between the Na atoms is maintained, ap- 
pearing however only in the upper half of the Na layer. 
The lower half is apparently cancelled by the Al surface 
states due to the anti-bonding coupling. A cross-section 
parallel to the surface passing through the electron den- 
sity maxima on top of the Na atoms is shown in Fig. 7c. 



We see also here a smeared out, free-electron- like distri- 
bution. 



3. Formation mechanism of the two surface state /resonance 

bands 

Combining the results above, we can conclude that the 
two bands are resulting from a coupling between the low- 
est lying (3s-derived) band of a free c(2 x 2) monolayer 
of Na and the surface state band of the clean Al (001) 
surface. Applying a simple two-term perturbation, thc.- 
ory, as commonly done in molecular-orbital theoryE3E3, 
this may be understood as being due to the formation 
of bonding and anti-bonding states, leading to the down- 
ward shift of the Al-derived band by 0.9 eV (calc.) with 
an increase in population and the upward shift of the Na- 
derived band by 0.8 eV (calc.) with a decrease in pop- 
ulation. Alternatively, it may be understood, as in the 
case of ionic crystals, that the shifts occur as the result 
of Coulomb fields between the two oppositely charged 
ionic layers. The ionization is to be expected due to the 
cationic nature of Na, which donates electronic charge to 
the aluminum. 
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FIG. 5. Free c(2 x 2)-Na layer: (a) band structure (b) 
charge density distribution of valence elctrons. The planes 
of cross-sections are denned in Fig. 4b. 

We show in Fig. 9 the total charge density, ptot., the 
density difference between the adsorption system and the 
corresponding Al surface (for which the positions of the 
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FIG. 6. Charge density distribution of the surface 
state/resonance of the main bands at V for (a) the clean sur- 
face, (b) the Al-derived band of the LT phase (hollow site), 
and (c) the Na-derived band of the LT phase. The positions 
of cross-sections are indicated in Figs. 4a and 4b. Large dots 
denote Na atoms and small dots Al atoms. The contours be- 
gin at 1.0 and the contour spacing is 1.0. The units are x 10~ 3 
e bohr -3 . 

Al atoms are kept at those of the adsorption system), 
Apdifi. = Ptot. — Pa\, an d the redistribution of charge, 
Ap rc dis. = Ptot. - pAi - P c (2x2)-Na, which subtracts out 
the charge density of the free Na monolayer (Fig. 5b), 
showing exactly where charge has been enhanced and 
depleted. It can clearly be seen from Apdiff. and A/? re di s . 
that charge enhancement occurs primarily at the maxima 
of the surface states, and indeed almost proportionally. 
We can also note in Ap re di s . some regions of depletion, 
showing that the electron charge has been transferred 
from the upper half of the region between the Na atoms, 
where the density for the free Na layer (see Fig. 5b) 
is much larger than that of the Al surface states (see 
Fig. 6a) . We may thus conclude that the electron transfer 
is occurring from Na atoms directly into the pre-existing 
surface states of Al. 

It is noted in passing that this character of the LT 
phase showg_.a_close analogy to the case of c(2 x 2)- 
Cs/W(001)BB, both having the c(2 x 2) four-fold hol- 
low structure. Although the electronic structure is much 
more complicated for Cs/W, the essential feature of the 
charge transfer is the same. Thus, the maxima of the 
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FIG. 7. Charge density distribution of the main surface 
state/resonances in a plane parallel to the surface passing 
through the electron density maxima (not the atom centres): 
The clean surface (a), the Al-derived (b) and Na-derived (c) 
bands of the LT (hollow) structure, respectively, and the va- 
cancy structure (d), the Al-derived (e) and Na-derived (f) 
bands of the RT (substitutional) structu re, respectively. The 
dots represent the positions of the Na atoms. The contours 
begin at 1.0 and the contour spacing is 1.0. The units are 
xl0~ 3 e bohr -3 . 

4d-derived surface states of the clean W(001) surface 
lie on top of the surface W atoms (see Fig. 7 of Ref.EJ). 
The maxima maintain their form upon Cs adsorption. 
Charge transfer takes place from Cs to these surface 
states. Hence, this case and ours may be regarded as 
representatives of the alkali-metal on metal systems for 
which the surface states of the substrate play an essential 
role. 

The upper Na-derived band is crossing the Fermi level 
and partly occupied, remaining "metallic" . In Figs. 6c 
and 7c its free-electron-like character can be seen. We 
see in Figs. 6b and 6c that the charge density of the Al- 
derived band lies well below the smeared out density of 
the Na-derived band, and also below the position of the 
Na atoms. Thus, the traditional picture of a thin metal- 
lic film covering an originally metallic substrate remains 
qualitatively valid. We see below that this is not the case 
for the RT phase. 

It is to be noted that the apparently weak influence of 
the Na adsorption on the Al surface states is limited to 
the high coverage of 0=0.5 for the c(2 x 2) structure, for 
which the density maxima for the Na valence electrons 
(Fig. 5b) lie between the Na atoms and situated just at 
the same site as the surface-state maxima. In fact, it has 
been found experimentallyliJ that at low coverages O=0 
to 0.15 the surface states are deteriorated by Na adsorp- 
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tion. Also it has been shown by DFT calculations for 
a Active ordered-structure model, p{2 x 2) with coverage 
0=1/4, that the surface-state maxima are moved from 
the on-top sites to the bonding-line positions between the 
Na and Al atoms. 
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FIG. 8. Same as Fig. 6 for the RT phase (substitutional 
site). The positions of cross-sections are specified in Figs. 4c 
and 4d. 



C. RT c(2 x 2) substitutional structure 

1. Al-derived band 

Figure 2c shows the experimental (left panel) and cal- 
culated (right panel) surface state/resonance bands for 
the RT c(2 x 2) substitutional structure. It can be 
seen that the experimental surface state/resonance band 
clearly exhibits a c(2 x 2) periodicity, in contrast to the 
LT (hollow) structure. As we will see below, this is due 
to the significantly reconstructed Al (001) surface. The 
calculations (right panel) show good agreement with ex- 
periment and indicate that the band is mainly Al-derived 
(circles) but has a small participation of Na (squares). 
The energy position of the main theoretical and experi- 
mental bands at T are 2.14 and 2.31 eV, respectively. 

The experimental results display some asymmetry in 
the intensity near k|j=0.5 in T-M, i.e. a weakening on 



the left side and a strengthening and upward shift on the 
right. A few alternative explanations may be possible 
for this anomaly. Among others, it may be related to 
the closing of the bulk band gap in this region, causing 
the state to go from being a pure surface state to a sur- 
face resonance (see Fig. 10c of Refill). Also, additional 
coupling to bulk states may become possible by surface 
Umklapp processes. The peak profile may become broad 
and complicated, as indicated by the split structure of 
the band near kn =0.5 in T-M. Different profiles for the 
left- and right-half may result from this broadened peak 
in the formation of the matrix element given by Eq. 4, 
the final state being different. 

An important finding in the present analysis is that 
the main band in the RT structure can be regarded as 
being derived from the surface state/resonance band of 
a Active, reconstructed clean Al surface, that is, the "va- 
cancy" structure, see Fig. 4c. In this structure the Na 
atoms are replaced by "vacancies" of the Al atoms. The 
surface state/resonance band structure of the vacancy 
structure is displayed in Fig. 3c. By comparison with 
Fig. 3b it is clearly seen that the main Al-derived band of 
the RT phase originates from that of the vacancy struc- 
ture and is only somewhat shifted down in energy (by 
0.81 eV) due to Na adsorption. 

Figures 8a and 8b show the single-state charge distri- 
bution at r for the main bands of the "vacancy struc- 
ture" and the RT phase, respectively. Figures 7d and 7e 
show their cross-sections parallel to the surface passing 
through the electron density maxima. Similarly to the 
clean Al(001) surface, the maxima of charge density for 
the vacancy structure lie on top of the uppermost sur- 
face Al atoms, having this time the c(2 x 2) periodicity. 
For the Al-derived band of the RT phase, the maxima 
also lie approximately at the same position as that of the 
vacancy structure with relatively small changes in their 
form. It can be noted that these states are rather strongly 
localized. This explains the well developed c(2 x 2) char- 
acter of the band found in Fig. 2c, and its relatively small 
dispersion (larger value of m*, see Table II). 

Similarly to the LT phase, we find in Fig. 2c a few 
additional theoretical bulk bands which are not present 
for the clean surface. These are apparently introduced 
by the c(2 x 2) periodicity as surface Umklapp processes. 
There is also indication of a third, relatively weak band 
lying at Ep and around T and M in both the experimental 
and calculated band structures. From our analysis (not 
shown), this band can be regarded as being derived from 
a second, relatively weak, surface state/resonance band 
of the vacancy structure. 

2. N a- derived band 

In Fig. 3b we see that the Na-derived band (filled 
squares) is found only in the calculation, lying completely 
in the positive energy range, about 2.7 eV higher in en- 
ergy than for the free c(2 x 2)-Na layer. Notably, this shift 
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is much larger than the downward shift of the Al-derived 
band ( 0.8 eV). This will be discussed below. 

The single-state charge distribution of the unoccupied 
Na-derived state at T is shown in Fig. 8c, indicating again 
an antibonding character. Figures 8c and 7f show its 
smeared out, free-electrondike structure with the max- 
ima residing this time above the Na atoms. 




FIG. 9. Total charge density distribution, ptot., (left panel), 
density difference, Apdiff., (middle panel), and density redis- 
tribution, Aprodis., (right panel), of the LT phase (hollow site). 
The cross-section is in the (OOI)Na plane (see Fig. 4b). The 
units are 10~ 3 e bohr~ 3 . In the left panel the first contour be- 
gins at 4.0 with a spacing of 4.0, for the middle panel the first 
contour begins at 0.6 with a spacing of 0.6, and in the right 
panel the contours are the same as the middle panel with the 
addition of a negative contour line (unshaded) at —0.6. 



3. Formation mechanism of the two bands 

We conclude that the two main bands are resulting, 
just as in the case of the LT phase, from the coupling 
between the states of a free c(2 x 2)-Na monolayer and 
the surface state/resonance of the Al surface, i.e., the 
vacancy structure in this case. 

The electron transfer can also be regarded as occur- 
ring from Na atoms into the surface state/resonance of 
the vacancy structure. This can be seen clearly in com- 
paring Fig. 10 with Fig. 8a. In the charge redistribution 
Apredis. in Fig. 10 we note some regions of depletion: the 
electron charge has been transferred in this case mainly 
from the region on top of the Na atoms, where the den- 
sity of the free Na layer dominates the density of the 
surface state/resonance of the vacancy structure. The 
charge depletion found in Fig. 10 also near the centres of 
the Al atoms may be interpreted as a result of an upward 
shift of the maximum of the surface states induced by Na 
adsorption. 



Aot. 


Pdiff. 


Presid. 















FIG. 10. Same as Fig. 9 for the RT phase (substitutional 
site). For the middle panel there is one contour line at —1.2, 
and in the right panel there are three negative contour lines 
(unshaded) at —1.8, —1.2, and —0.6. 

These findings, combined with the unbalancedly large 
upward shift of the Na-derived band and the existence 
of the third peak at Ep, both mentioned above, indicate 
that the coupling between the free Na monolayer states 
and the Al surface states/resonances takes place in a lit- 
tle more complicated way than in the LT case. It can no 
longer be interpreted in a simple two-term-coupling per- 
turbation scheme. Apparently the coupling is too strong, 
and involves other states, namely, the p z states of Na ly- 
ing originally at +1 eV (Fig. 5a) and the second surface 
states/resonances of Al lying originally at eV (Fig. 3c). 

We have seen above that for the LT phase the tradi- 
tional picture of a thin metallic film covering a metallic 
substrate remains qualitatively valid. For the RT phase, 
on the other hand, we see in Fig. 8b that the completely 
filled Al-derived band has the maxima lying at the same 
height as the Na atoms. The maxima of the Na-derived 
band lie indeed higher, but they are empty. Thus, as 
far as the electronic states in the energy range near Ep 
are concerned, the Al-derived band constitutes practi- 
cally the "surface" of this structure. The Na monolayer 
(or the Na/Al composite monolayer) cannot be regarded 
as a metallic film on a metallic substrate; rather, ex- 
cept for the still existent background bulk-continuum, 
the electronic structure of the monolayer may be viewed 
as being analogous to that of an ionic crystal, like NaCl, 
where the completely filled Al-derived band replaces the 
valence band of CI, and forms with the empty Na-derived 
band a "band-gap" of about 3.5 eV across the Fermi en- 
ergy. 

We note in passing that we find a similarity of the 
surface electronic structure of the RT phase to that of 
the system (2 x 2)-2Na/Al (111) which forms a composite 
double layer surface alloy with a similar, but more £QgL- 
plex intermixing of Na and Al in the surface layersOO. 
For this system we find also a filled Al-derived band orig- 
inating from the corresponding "vacancy structure" , and 
an empty Na-derived band. (The-Jatter may be assigned 
to that found by Heskett et al. using inverse photoe- 
mission spectroscopy.) The similarity indicates that also 
this layer cannot be regarded simply as a thin metallic 
layer. 



V. DISCUSSION 

For all three systems studied (Al(001) and the two 
c(2 x 2)-Na/Al(001) phases), we find good over-all agree- 
ment of the surface state/resonance band structure be- 
tween experiment and theory. In each case there is 
a prominent Al-derived surface-state/resonance band, 
showing similar free-electron-like parabolic dispersion at 
the band bottom near T, as indicated by broken lines in 
Fig. 2. We find that the observed and calculated val- 
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ues of the Fermi wave vector, kp, and the effective mass, 
to* of the main band (Tab. II) agree only roughly, prob- 
ably due to the approximations used in calculation, as 
already mentioned. In any case these quantities are to 
be regarded as global parameters specifying only the ge- 
ometry of the bands, and are not intended to indicate 
that the bands are free-electron-like. In fact, the wave 
functions of the Al-derived bands are rather strongly lo- 
calized, as we see from the charge density distributions. 
The dispersion of the bands may be regarded as the result 
of the overlap between the localized wave functions. 

The picture of charge transfer taking place from Na to 
the surface state/resonances of Al is also supported by 
the almost equal values of the work function change.-(d|6-| 
crease) A$ w-1.6 eV, obtained both experimentallyE§EZl 
and theoretically in the present work for the LT and RT 
phases. Obviously the value of A$ results from p re dis.) 
shown in Figs. 9 and 10, as the z-component of the dipole 
moment. As already mentioned, the minimum and max- 
ima of p r odis. correspond to the maxima of the Na- and 
Al-derived states, respectively. This can be seen by com- 
parison of Figs. 9 and 10 with Figs. 6b, c and Figs. 8b, c. In 
these figures, the vertical distances between the maxima 
of the Na- and Al-derived states are found to have the 
ratio of about 1.0 to 0.6 between the LT and RT phases. 
On the other hand, a comparison between Figs. 5a and 
2b, c shows that the decrease of the occupancy of the Na- 
derived states amounts to ca. 0.6 for the LT phase, in 
contrast to 1.0 for the RT phase (being empty). This 
lead to the same amount of A<f> between the two phases, 
as a product of charge and distance. It is essential in this 
consideration to note the fact that, while the forms of the 
Al-derived states remain always almost unchanged, the 
forms of the Na-derived states change drastically, from 
that of Figs. 5b to those of Figs. 6c and 8c. 

Our results verify the formation of the two surface 
state/resonance bands. Various experimental methods 
of studying surface electronic properties other than ul- 
traviolet photoelecton spectroscopy, as used here, may 
be useful for finding out characteristic features induced 
by these two bandp_.£articularly, EELS (electron energy 
loss spectroscooy-J)l3lEl and IPS (inverse photoemission 
spectroscopy )l3cj may be useful for studying the effects 
of the bands lying in the energy range above the Fermi 
level. 

For the RT phase, the fact that the surface band struc- 
ture is constituted of a filled band and an empty band 
with a gap of approximately 3.5 eV, may play an im- 
portant role in various properties of the surface. We 
may think of, in addition to EELS, an anomalous fea- 
ture in optical reflection spectracJ near 3.5 eV and a cor- 
responding, anomaly in the dispersion curve of surface 
plasmonsEj. 

A significant effect may be expected in various 
surface-sensitive methods,— such as ion neutralization^}, 
metastable deexcitationOEj etc., which would reflect the 
dominance of the surface states/resonances in the out- 
ermost surface region, and hence expose the difference 
in the character of these states between the LT and RT 



phases. We may also think of the relevance of the occp-i 
pancy of surface statesYresonances on surface diffusions 
and catalytic activityEl 

VI. CONCLUSION 

We have analysed the electronic structure of the 
metastable hollow and stable substitutional structures 
of c(2 x 2)-Na on Al(001). The calculated surface 
state/resonance bands agree well with those measured 
by angle-resolved photoemission experiments. It is found 
that in both phases, two pronounced bands appear as the 
result of a characteristic coupling between the valence- 
state bands of a free c(2 x 2)-Na layer and the surface- 
state/resonance bands of the corresponding (i.e. clean 
and reconstructed) Al surfaces. While the experimen- 
tal band structure of the substitutional structure shows 
a clear c(2 x 2) character due to the significant recon- 
struction of the surface Al layer, the hollow structure 
does not: The main band, in fact, exhibits a quasi- 
(1 x 1) periodicity like that of the clean surface. For the 
stable substitutional structure, the unoccupied surface 
state/resonance band lies completely above the Fermi 
level, leading to the formation of a surface-state band 
structure that resembles that of an ionic crystal. We 
await experimental verification of the predicted unoccu- 
pied surface state/resonance band, and of the difference 
in the properties of the LT and RT phases in relation to 
the character of the surface states/resonances. 
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vacancy 


big circles 
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small circles 


Al 


0.3 
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TABLE I. Legends of the 


symbols used in calculated 2D 


band structures in Figs. 2 and 3. The "minimum fraction" is 


defined in 


the text (see HI Calculation method) . 
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